Scaffold-mediated Nucleation of Protein Signaling Complexes: Elementary Principles 
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Proteins with multiple binding sites play important roles in cell signaling systems by nucleating 
protein complexes in which, for example, enzymes and substrates are co-localized. Proteins that 
specialize in this function are called by a variety names, including adapter, linker and scaffold. 
Scaffold-mediated nucleation of protein complexes can be either constitutive or induced. Induced 
nucleation is commonly mediated by a docking site on a scaffold that is activated by phosphoryla- 
tion. Ifere, by considering minimalist mathematical models, which recapitulate scaffold effects seen 
in more mechanistically detailed models, we obtain analytical and numerical results that provide 
insights into scaffold function. These results elucidate how recruitment of a pair of ligands to a 
scaffold depends on the concentrations of the ligands, on the binding constants for ligand-scaffold 
interactions, on binding cooperativity, and on the milieu of the scaffold, as ligand recruitment is 
affected by competitive ligands and decoy receptors. For the case of a bivalent scaffold, we obtain 
an expression for the unique scaffold concentration that maximally recruits a pair of monovalent 
ligands. Through simulations, we demonstrate that a bivalent scaffold can nucleate distinct sets 
of ligands to equivalent extents when the scaffold is present at different concentrations. Thus, the 
function of a scaffold can potentially change qualitatively with a change in copy number. We also 
demonstrate how a scaffold can change the catalytic efflciency of an enzyme and the sensitivity 
of the rate of reaction to substrate concentration. The results presented here should be useful for 
understanding scaffold function and for engineering scaffolds to have desired properties. 
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effect, combinatorial inhibition, adaptor 



I. INTRODUCTION 



Proteins involved in cell signaling typically possess 
multiple binding sites, and as a result, a common fea- 
ture of these proteins is their ability to interact with sev- 
eral binding partners to form heterogeneous protein com- 
plexes pp. Some proteins specialize in the nucleation of 
protein complexes [2]. For example, Ste5 is a well-studied 
scaffold protein that acts in a yeast MAP kinase cascade; 
it binds all three kinases of the cascade, Stell, Ste7, and 
Fus3 [3 . The interaction of FusS with Ste5 is modulated 
by multisite phosphorylation of Ste5 [3] . One function of 
Ste5 is to co-localize enzyme-substrate pairs (e.g., Stell 
and Ste7) in the cascade [5] , which contributes to switch- 
like cellular responses to signals 4J. Two well-known 
examples of scaffold-like proteins that nucleate protein 
complexes in mammalian signaling systems are the Grb2 
adapter protein in the epidermal growth factor receptor 
(EGFR) pathway [Ij and the LAT linker prote in in the 
T cell receptor (TCR) pathway [7HS]. The binding sites 
in LAT are activated by tyrosine phosphorylation, and 
as a result, the valence of LAT and the abundance of 
functional LAT binding sites can be quickly changed by 
changes in kinase and phosphatase activities, with these 
changes serving to modulate the effect of LAT on signal- 
ing [To]. In addition to proteins with dedicated scaffold 
functions, proteins that have catalytic activities and pro- 



tein complexes, such as ligand-induced receptor dimers, 
can also function as scaffolds, i.e., take on the function of 
facilitating the formation of multicomponent complexes 
through multivalent binding. 

Many functions have been ascribed to scaffold pro- 
teins [21 [9l [TTUI4] . For example, scaffolds have been 
suggested to control the specificities of enzymes, coor- 
dinate spatiotemporal aspects of signaling, and amplify 
or attenuate signals. To better understand the mech- 
anisms by which scaffolds influence cell signaling, re- 
searchers have formulated mathematical models to study 
the effects of scaffolds on specific signaling systems and 
to study the generic properties of scaffolds [H fT51l22) . 
The ability of a scaffold to either amplify or attenuate 
signaling was demonstrated by Levchenko et al. [TB] us- 
ing a mathematical model for a MAP kinase cascade. 
In this model, scaffold enhancement of signaling results 
from scaffold-mediated nucleation of an enzyme with 
its substrate, whereas scaffold inhibition of signaling re- 
sults from excess scaffold, which separates enzyme and 
substrate into distinct enzyme-scaffold and substrate- 
scaffold complexes. This effect is equivalent to the well- 
known prozone effect in antibody-antigen reactions [15j . 
Heinrich ct al. ,17J further studied these effects of scaf- 
folds on protein kinase cascades. More recently, Locasale 
et al. [20 studied how scaffolds affect competition be- 
tween phosphatases and kinases, finding that scaffolds 
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can have positive and negative effects on signaling for 
reasons related to the kinetics of signaling events. 

Modeling has improved our understanding of scaffold 
function, but many of the reported modeling studies of 
scaffold function have been based solely on simulations, 
which makes the results dependent on the parameter val- 
ues considered in simulations. In addition, many of the 
models used to study scaffolds incorporate mechanistic 
details that are likely to be extraneous for understand- 
ing scaffold function broadly. Greater theoretical under- 
standing of the design principles of scaffolds could per- 
haps be obtained by focusing on minimalist models that 
are relatively easy to analyze. With a better understand- 
ing of scaffold design principles, we can hope to manip- 
ulate scaffold function through precise tuning of scaffold 
and system properties. Manipulating the properties of a 
scaffold or its milieu to alter the behavior of a cell sig- 
naling system is known to be feasible [S] ESI [21] • For ex- 
ample, Lim and co-workers have demonstrated that the 
response of a cell to a signal can be changed qualitatively 
by altering the binding specificities of a scaffold [5]. 

Here, by considering minimalist models for multiva- 
lent scaffold-ligand interaction, we investigate the re- 
cruitment of ligands to a scaffold and the effects of a 
scaffold on an enzymatic reaction. We obtain a number 
of analytical and numerical results, mostly for a model 
that characterizes the equilibrium formation of a ternary 
complex composed of a scaffold and two binding part- 
ners. In this ternary-complex model, the scaffold is biva- 
lent and and its monovalent binding partners are taken 
to be an enzyme-substrate pair. Each binding partner, 
or ligand, interacts with one of the two binding sites of 
the scaffold. We evaluate the relevance of the ternary- 
complex model through comparisons with more mecha- 
nistically detailed models for cell signaling systems that 
involve scaffolds. The ternary-complex model is similar 
in mathematical form to a variety of models that have 
been used to study receptor signaling [521 HSI a-nd multi- 
valent ligand- receptor binding [27-25]. 

The remainder of this report is organized as follows. 
In Section [ll| we introduce models and present analytical 
results relevant for understanding the effects of scaffold 
valence, negative and positive cooperativity in binding 
of multiple ligands to a scaffold, competitive inhibitors of 
ligand-scaffold interactions, and decoy receptors. We also 
develop an approximate rate law for scaffold-mediated 
enzyme kinetics. Our most significant analytical result 
is a design equation that characterizes how the scaffold 
concentration that maximally nucleates a ternary com- 
plex depends on equilibrium binding constants and lig- 
and concentrations. Interestingly, this scaffold concen- 
tration is independent of cooperativity in ligand binding. 
In Section [III[ we present numerical results that comple- 
ment the analytical results. The numerical results in- 
clude a demonstration that our simple ternary-complex 
model recapitulates scaffold effects seen in more mech- 
anistically detailed models and a demonstration that a 
single scaffold can nucleate distinct signaling complexes 
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FIG. 1. Minimalist models for ligand-scafTold interactions. 
(A) A scaffold with n sites (5*1, . . . , Sn)- Each site Si interacts 
with a distinct ligand L,;. (B) Ternary-complex model. A 
bivalent scaffold with sites a and h interacts with monovalent 
ligands A and B, respectively. (C) A scaffold site a interacts 
with a ligand A and a competitive inhibitor /. The ligand 
A interacts with a decoy receptor D. (D) A reaction scheme 
for a scaffold-mediated enzymatic reaction that produces a 
product A! . The substrate is scaffold ligand A and the enzyme 
is scaffold ligand B. See text for discussion. 



at different scaffold concentrations. We conclude with a 



brief discussion (Section IV I . 



II. MODELS AND ANALYTICAL RESULTS 

The models that we will consider are illustrated in 
Fig. [T] and described below, along with various analyt- 
ical results relevant for understanding scaffold function. 



A. Multivalent scaffold with independent binding 
sites 

We consider a scaffold protein with n > 2 binding sites. 
Each scaffold binding site i interacts independently with 
a distinct monovalent ligand i, as illustrated in Fig.[T]^A). 
Thus, we consider the following reactions: 



Si ^ Li \ Bi^ i — 1, 



(1) 



where Si represents free site i, Li represents free ligand 
i, and Bi represents bound site i (or equivalently, bound 
ligand i). Using the same symbols for chemical species 
in Eq. ([I]) for the corresponding concentrations, we can 
write the following equilibrium relations, which are de- 
rived from the law of mass action: 



KiBi ~ SiLi, i — I, 



(2) 



where Ki is an equilibrium dissociation constant. As- 
suming conservation of mass, such that Sq = Si -\- Bi 
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and = Li + Bi, where Sq is the total concentration 
of scaffold protein and Li^ is the total concentration of 
ligand i, we find from Eq. ([2| that 

B^ -{So + L,„+K,)B, + SoL,„=0 . (3) 

From the quadratic formula and the physical constraint 
that Bi < min(iif, , 5*0), it further follows that 



B, = 



So + ^{So + L,^,+K,)^ -ASoL,, 



(4) 

This equation and the mass conservation relations given 
above can be used to completely determine the equilib- 
rium state of the system described in Fig. [I|A). 

We will now show that there is a unique scaffold con- 
centration S'gP* that maximizes C„, the equilibrium con- 
centration of scaffold protein with all n sites occupied. 
The analytical results that follow complement numerical 
results obtained in earlier modeling studies of scaffold 
function [nillMI]- 

We note that the equilibrium quantity Bi/So can be 
interpreted as the probability that scaffold binding site 
i is occupied at equilibrium. Likewise, C„/5o can be 
interpreted as the joint probability that all n scaffold 
binding sites are occupied. Thus, because each scaffold 
binding site is assumed to interact with its cognate ligand 
independently, we can write 



nl — n 

n — '^0 



(5) 



Trivially, for the case 5*0 = 0, C„ = 0. In other words, 
there is no complex in the absence of scaffold. Also, from 
Eq. ([5|, we can see that limS(,_>.oo Cn — 0, because each 
Bi is finite. In fact Bi can be no greater than Li^ . From 
inspection of Eqs. Q and ([s]), one can see that C„ is a 
continuous and differential function of So- Thus, from the 
mean value theorem, the derivative of C„ with respect to 
So must vanish at some finite value of 5*0, which we will 
denote as S^^^^ . In other words, 



dSo 



= 



(6) 



Sn—Sn 



Because Sq^*' is unique, as we will show below, and C„ 
is positive for all finite 5*0 (on physical grounds), C„ is 
maximum at 5*0 = Sq^^ . 

To see that Sq^^ is unique, consider the following ex- 
pression, which is derived from Eq. ^ by differentiating 
Cn with respect to 5*0: 



So dCn sr-^ , , X 

i—i 



where 



So dBi So{Bi — Lig 



Bi dSo B- — SoLi^ 



(7) 



(8) 



Note that g„ is the slope of the plot of log C„ vs. log 5*0. 
The expression for bi was obtained by differentiating each 
term in Eq. ([3| with respect to So and then solving for 
dBi/dSo- From inspection of Eq. ([s]), one can see that 
< 6i < 1 for finite So- In fact, it can be shown that 
lims^^o bi = 1 and lim5,-|_j.oo bi — 0. From these consid- 
erations and inspection of Eq. ([?]), one can see that Sq^^ 
is unique if bi decreases monotonically as Sq increases for 
all i. Differentiating the expression for bi given in Eq. ([S]) 
with respect to 5*0, we find 



dbi 
dSo 



Bi{Bi 



(Bf- 



~ So)iBf + SpLig) 



< 



(9) 



Note that the terms Bi — Li^ , Bi — So, and Bf — SoLi^ are 
all negative. Thus, each hi is a monotonically decreasing 



;opt 



function, as is gn by extension, and as a result, Sq 
unique. 

The results presented above provide insight into the 
shape of the log-log scaffold dose-response curve, mean- 
ing the plot of logC„ vs. logS'g. Recall that the slope 
of this curve, which is concave, is given by Eq. ^. As 
can be seen from Eq. ([t]), for Sq < Sq^^, 1 > 5„ > 0, 
whereas for Sq > Sq^^, > (;„ > 1 — n. For n = 2, 1 
approaches 1 as So approaches either or 00. However, 
for n > 3, the two asymptotic limits of are different, 
with the curve becoming steeper at some point in the 
excess scaffold regime (^o > S°^') than at any point in 
the regime where So < Sq^^- Thus, the log- log scaffold 
dose-response curve is asymmetric. 

To gain insight into the effect of scaffold valence on 
scaffold function, measured by nucleation of a complex 
in which all scaffold binding sites are occupied, let us 
consider two scaffolds in the same milieu with the same 
properties except that one scaffold has an extra bind- 
ing site n and cognate binding partner L„. We impose 
this equivalence constraint to ensure a fair comparison 
[30l I3T] . We will denote the valences of the two scaffolds 
as TO > 2 and n = to -I- 1. From Eq. ([5| and our equiv- 
alence constraint, it follows that C„ = CmBn/So and 
further that C„ < Cm for any given scaffold concentra- 
tion. Thus, with all factors being as equal as possible, 
scaffold-mediated nucleation of C„ is less efficient than 
nucleation of Cm- Another inherent effect of valence on 
scaffold function is greater sensitivity of C„ to changes 
in scaffold concentration in the excess scaffold regime. 
From Eq. ([7]), for So sufficiently large such that < 
and gn < 0, it can be seen that gn is larger in magnitude 
than gm- This result is consistent with our earlier discus- 
sion of asymmetry in the log-log scaffold dose-response 
curve. 



B. Bivalent scaffold with cooperative ligand 
binding 

We will now consider a bivalent scaffold that interacts 
with a pair of monovalent ligands A and B (Fig. [iJB)). 
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Ligand A interacts with a scaffold binding site a, and 
ligand B interacts with a scaffold binding site b. Equa- 
tion ([T]) with n — 2 applies, but we will not use the 
nomenclature of Eq. ([T]) in our treatment of this spe- 
cial case, which because of its simplicity will allow us to 
obtain an analytical expression for the scaffold concen- 
tration that maximizes formation of the ternary complex 
composed of A, B and the scaffold. Nucleation of this 
complex can be taken as a measure of the functional ac- 
tivity of the scaffold if A and B constitute an enzyme and 
substrate or if instead the two ligands constitute func- 
tionally related enzymes recruited to the scaffold to act 
in concert, as in substrate channeling [32] ■ We will re- 
lax the assumption that ligands interact with the scaffold 
independently and allow for cooperative binding. Thus, 
using the law of mass action, we write the following equi- 
librium relations: 



KaCa — SfAf , 

KbCb — SfBf , 

KaCab = 4>AfCb 



(10) 

(11) 
(12) 



where Sf is the concentration of free scaffold, Af and Bf 
are the concentrations of free ligands A and -B, Ka is an 
equilibrium dissociation constant that characterizes a-A 
interaction, Ca is the concentration of the binary com- 
plex composed of A and the scaffold, Kh is an equilibrium 
dissociation constant that characterizes b-B interaction, 
Cb is the concentration of the binary complex composed 
of B and the scaffold, and Cab is the concentration of 
the ternary complex composed of A, B, and the scaf- 
fold. The factor 0, which is positive and dimensionless, 
is introduced to characterize cooperativity in binding of 
A and B to the scaffold: </> < 1 indicates negative co- 
operativity, (/) > 1 indicates positive cooperativity, and 
(j) — 1 indicates that the ligands bind the scaffold inde- 
pendently. We will assume conservation of mass, such 
that the following relations hold: 



Sf) — Sf + Ca - 

Aa = Af + Ca 
Bo = B 



f 



C, 



Cb + Cab 
' Cab J 
Cab ■ 



(13) 
(14) 
(15) 



where Sq is the total scaffold concentration, Aq is the 
total concentration of ligand A, and Bq is the total con- 
centration of ligand B. 

We will now derive an expression for Sq^*" , the 
total scaffold concentration that maximizes Cab, or 
equivalently, the total scaffold concentration for which 
dCab/dSo = 0. From Eqs. pOl) and (Il2l, we find 



Cab 



CaCb 



(16) 



We note that, for 0=1, Cab = CaCb/Sf = {Ca + 
Cab){Cb + Cab)/So, which is equivalent to Eq. ^ with 
n = 2. From Eqs. (10) and (14), we find 



Ca 



Sfi^O ~ Cab) 
Sf+Ka 



(17) 



Similarly, from Eqs. (11) and (15), we find 



Cb 



SfiBo — Cab) 

Sf + Kb 



From Eqs. (16)H18), we find 



Cab ~ RCab + AqBq = 



(18) 



(19) 



(20) 



Differentiating each term in Eq. ( |19[ ) with respect to 5*0, 
we find 



where 



R=Ao + Bo + {Sf + Ka){Sf + Kb)/{cj>Sf) 



dCab _ f Cab (t) 
dSo ~ UCa. 




(21) 



From this equation, one can see that the derivative 
dCab/dSo uniquely vanishes at Sf = ^/KaKb- Thus, S?^ 
is the value of 5*0 at which Sf = \/KaKb- Using Eqs. ( 17 ) 
and (18), we can rewrite Eq. (13) as follows: 



5o — Sf + - 



SfA 



f^O 



SfBo 



1 - 



S 



f 



S 



f 



Sf +Ka Sf+ Kb 

(22) 



Sf+Ka Sf+Kb 

Substituting y/KJ<b for Sf in Eq. (22), we find 

^0 



= ./KJ{~b 



Bn 



l + ^/KjK~b l+./K^a 



(23) 



Notably, S'^^'^ is independent of the cooperativity factor 
(f). For the special case Ka — Kb — Kb, Ku + 

{Aq + Bq) /2. For the special case Aq = B() = Lq, S^^*^ — 

y/KaKb + Lo. 

It is straightforward to find an expression for the max- 
imum concentration of scaffold-nucleated ternary com- 
plex, which we will denote as CJ^^. We simply use the 



quadratic formula to solve Eq. ( 19 ) for the case where 

This expression for 



R = Aq + Bp + (j)-^{VK; + V K^ 
R is obtained by substituting ^/KaKb for Sf in Eq. (20) 
Unlike ST, CI 



opt ^max (igpgu(jg on (j) (bccause R depends on 
For the special case where (/> = 1, Ka = Kb = Kn and 



Ao - So = Lo, we find = WKd + io 



C. Effect of a competitor or decoy receptor 

The in vivo milieu of a scaffold is complex. A given 
binding site on the scaffold may interact with a specific 
ligand as well as a number of competitors. Likewise, the 
ligand may interact a number of binding partners be- 
sides the scaffold. To study the effects of the milieu on 
the function of a scaffold, we consider the simple scenario 
illustrated in Fig. [l](C). A ligand A interacts with a scaf- 
fold binding site a in competition with an inhibitor /. 



Cab 
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The hgand also interacts with a decoy receptor D. Thus, 
we consider the following equilibrium relations: 

KaBa ^SaAf , (24) 

K,B, = Salf , (25) 
KdBd = DfAf, (26) 

where Sa is the concentration of free scaffold binding site 
a, Af is the concentration of free ligand A, Ba '\& the con- 
centration of scaffold bound to ligand A, // is the con- 
centration of free inhibitor /, Be is the concentration of 
scaffold bound to competitive inhibitor I, Df is the con- 
centration of free decoy receptor D, and Bd is the concen- 
tration of ligand A bound to decoy receptor D. The equi- 
librium dissociation constants Ka , and Kd character- 
ize a- A interaction, a-I interaction and A-D interaction, 
respectively. We will assume that the following mass con- 
servation relations hold: Sq = Sa-\- Ba + Be, where is 
the total scaffold concentration, and Aq — Af + Ba + Bd, 
where Aq is the total ligand concentration. 

To quantify the effect of a competitive inhibitor or de- 
coy receptor on interaction between a scaffold binding 
site and a ligand, we can find the apparent equilibrium 
dissociation constant K'^ satisfying the relation K'^^Ba = 
S'^A'^ (cf. Eq. g), where S', = So - Ba = Sa + B, is 
the apparent concentration of free scaffold binding site a 
and A'^ = Aq — Ba = Aj + Bd is the apparent concen- 
tration of free ligand A. From these considerations and 
Eqs. (p4l)-((26l), we find 



KL =Ka I 



II 

Ke 



1 



Kd 



(27) 



Thus, a competitive inhibitor or decoy receptor simply 
has the effect of reducing the apparent affinity of a ligand 
for its scaffold binding site. 

In the previous section, we provided equations that elu- 
cidate how S'^^ ^ and C™;^^ depend on binding constants. 
Given Eq. (27), we can now interpret these equations as 



also defining how competitive inhibitors and decoy recep- 
tors in the milieu of a scaffold affect S°^'^ and C™j^^. A 
competitive inhibitor or decoy receptor can shift the total 
scaffold concentration that maximally nucleates a signal- 
ing complex. For example, under certain conditions, a 
competitive inhibitor can promote scaffold-mediated nu- 
cleation of a signaling complex by reversing the effect of 
excess scaffold (i.e., the effect of 5*0 > Sq^*"). Thus, the 
function of a scaffold might be controlled by regulating 
the expression levels of competitive inhibitors and decoy 
receptors. It should also be noted that competitive in- 
hibitors and decoy receptors can obscure the relationship 
between S'^^'^ in vivo and binding parameters determined 
in vitro. 



D. Effect of a scaffold on the rate of an 
enzyme-catalyzed reaction 

To study the effect of a scaffold on the rate of an 
enzyme-catalyzed reaction, we consider the scheme pre- 



sented in Fig. [TJD). This scheme represents a sim- 
ple extension of the classical Michaelis-Menten reaction 
scheme: 



B ^ {AB) A' + B , 
k-i 



(28) 



where ^4 is a substrate, B is an enzyme, (AB) is an 
enzyme-substrate complex. A' is a product, and fci, 
and fccat are rate constants. As is well known, a 
pseudo steady-state assumption leads to the following 
approximate expression for the rate of product formation 
V=k,,tiAB): 



V ^ 



kcat BoAf 



Km+A 



f 



(29) 



where Af is the concentration of free substrate, Bq is the 
total concentration of enzyme, and Km = (fc_i-|-fccat)/fci- 
For simplicity, we will assume that fccat ^ fc-i such that 
Km ~ fc_i/fci. When Km is set equal to k^i/ki, it can 



be said that Eq. ( 29 ) depends on a rapid equilibrium 
assumption. 



To obtain an expression analogous to that of Eq. ( 29 ) 
for the reaction scheme of Fig. [iJD), we will assume 
that V , the rate of product formation when an en- 
zyme and substrate are co-localized on a scaffold, is 
proportional the concentration of the ternary complex 
of the scaffold, enzyme and substrate, which we will 
denote as Cat- In other words, we will assume that 
V = hkcs^tCab = hkcatiBo - Bf - Cb), where h > 
is a constant introduced to characterize the effect on re- 
activity of scaffold-controlled positioning of the enzyme 
and substrate. Positioning/orientation of an enzyme- 
substrate pair is a known function of some scaffolds [14j . 
Note that we have assumed conservation of enzyme, such 
that Bq = B f + Cb + Cab, where Bq is the total concen- 
tration of enzyme, Bf is the concentration of free en- 
zyme, and Cb is the concentration of the binary complex 
of enzyme and scaffold (Fig. [T]JD)). Under a rapid equi- 
librium assumption, we take the following relations to 
hold: KfjCb — SfBf and KaCab = AfCb, where Sf is the 
concentration of free scaffold and Ka and K^ are equi- 
librium dissociation constants. Note that Ka = k-a/ka 
and Kf, — k^f,/kb, where ka, k-a, k^ and k^i, are rate 
constants in the reaction scheme of Fig. [iJD). Thus, 



V 



hkca.tBoAf 



KL 



A 



f 



(30) 



where K'^ = Ka{l + K^/Sf). If h in Eq. ([30]) is ap- 
proximately 1 (i.e., scaffold effects on the reactivity of 
scaffold-tethered enzyme and substrate are inconsequen- 
tial, which can be the case ^ ) , then the effect of a scaffold 
on the rate of an enzyme-catalyzed reaction lies in the dif- 
ference between Km and K'm, the apparent Micheaelis- 
Menten constant for the scaffold-mediated reaction. Note 
that K'm is unrelated to Km and is determined by the 
affinities of the enzyme and substrate for binding sites 
on the scaffold and by the concentration of free scaffold. 
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Fraction of receptors in dimers 

FIG. 2. Comparison of predictions of the ternary-complex 
model and more mechanistically detailed models capturing 
the effects of scaffolds. (A) Comparison between the results of 
the ternary-complex model (solid line) and a model for a MAP 
kinase cascade incorporating a scaffold [16 . The broken line 
corresponds to the concentration of phosphorylated MAPK 
in Ref. [16]. Parameter values correspond to those used in 
Ref. [H]: Ao = 0.2 ^M, Bq = 0.4 and Ka = Kt = 
5 nM. (B) Comparison between the results of the ternary- 
complex model (solid line) and a model for early events in 
IgE receptor signaling |33l 134] . The broken line corresponds 
to the amount of Syk autophosphorylation predicted by the 
model of Ref. [25 for the case of 4 x 10^ Syk molecules per 
cell, as in Fig. 3 of Ref. [l] . Other parameter values are as in 
Ref. [34]. Note that the Kd for Syk binding to phosphorylated 
receptor is 2167 molecules per cell. 



Thus, a scaffold can profoundly affect the catalytic ef- 
ficiency of an enzyme even if the stereochemical effects 
of the scaffold on reactivity (captured in h) are inconse- 
quential. 



III. NUMERICAL RESULTS 
A. Relevance of simple ternary-complex model 

To investigate scaffold-mediated nucleation of signal- 
ing complexes, we will focus, except as noted, on a biva- 
lent scaffold that recruits a pair of monovalent ligands, 
i.e., on our ternary-complex model given by Eqs. |10[|15| 
and illustrated in Fig. [i|B). To demonstrate the rele- 
vance of this minimalist model, we will consider two more 



mechanistically detailed models for cell signaling systems 
reported in the literature. The first of these models char- 
acterizes the effect of a scaffold on a MAP kinase cas- 
cade [TB]. The second model characterizes early events 
of immunoreceptor signaling, in which a ligand-induced 
receptor dimcr nucleates a signaling complex [SS] [34] . 

In the model of Levchenko et al. fT^ , MAPK is phos- 
phorylated as a result of a cascade of binding and phos- 
phorylation reactions in which a scaffold co-localizes 
MAPK with its kinase, MAPKK. In Fig. [^JA, we plot 
the steady-state level of phosphorylation of MAPK pre- 
dicted by the model of Levchenko et al. [16] as a function 
of scaffold concentration. In addition, we plot the de- 
pendence of scaffold-nucleated ternary complex concen- 
tration (Cab) on scaffold concentration using the param- 
eters that characterize the binding reactions of MAPK, 
MAPKK and scaffold in the original model. As can be 
seen, if we assume that MAPK phosphorylation is re- 
lated to the amount of MAPK and MAPKK co-localized 
by scaffold, then our simple model closely predicts the 
steady-state level of MAPK phosphorylation predicted 
by the model of Levchenko et al. [16 even though the 
simple model omits many of the reactions included in 
the more mechanistically detailed model. This result 
suggests that the steady-state output of the MAP ki- 
nase cascade is dominated by co-localization of MAPK 
and MAPKK by the scaffold. We note that the MAPK- 
scaffold and MAPKK-scaffold equilibrium dissociation 
constants in the model of [IS] are each 5 nM. As a re- 
sult, Sq^^ is largely determined by the concentrations of 
MAPK (0.4 ^M) and MAPKK (0.2 ^M), as can be seen 
from Eq. ([23]). 

In the mechanistic model for early events in IgE recep- 
tor (FceRI) signaling developed by Goldstein et al. [S3] 
and Faeder et al. [34], a bivalent ligand induces recep- 
tor dimerization; receptor phosphorylation; recruitment 
of Syk, a cytosolic protein tyrosine kinase, to phospho- 
rylated receptors; and autophosphorylation of receptor- 
associated Syk. In this model, a phosphorylated recep- 
tor dimer acts as a scaffold that co-localizes two copies 
of Syk, enabling one copy of Syk to phosphorylate its 
neighboring copy of Syk. In Fig. [2^3, we plot the steady- 
state level of Syk autophosphorylation predicted by the 
detailed model with the parameters of Ref. [31] as a 
function of receptor dimer concentration. We also plot 
the dependence of Cab on phosphorylated receptor dimer 
(the effective scaffold) concentration using the parame- 
ters of Ref. [34 that characterize the binding reactions of 
phosphorylated receptor and Syk. As can be seen, if we 
assume that Syk autophosphorylation is related to the 
number of receptor dimers associated with two copies of 
Syk, then the ternary-complex model accurately predicts 
the amount of ligand-induced receptor dimers that yields 
maximum Syk autophosphorylation. This result further 
demonstrates that the minimalist ternary-complex model 
can provide insights into the effects of a scaffold on cell 
signaling. 

We caution that the ternary-complex model is not 
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guaranteed to recapitulate the system properties of more 
sophisticated, more mechanistic models. Our point is 
only that the model, despite its simplicity, can be use- 
ful for gaining insights into signaling events in which 
scaffold-mediated nucleation plays a role, as illustrated 
by the qualitative similarities of the two sets of curves 
in Fig. [2] The ternary-complex model can fail to cap- 
ture the effects of a scaffold in the context of a specific 
cell signaling system for myriad reasons. For example, 
an important caveat of the model is that it is an equi- 
librium model. Thus, if interactions of binding partners 
with a scaffold are slow relative to the time scale of pro- 
cesses affected by scaffold-mediated nucleation, then the 
ternary-complex model will fail to capture the true ef- 
fects of the scaffold, which will be dominated by the ki- 
netics of ligand-scaffold interactions. Also, the scope of 
the ternary-complex model is rather circumspect. For 
example, it does not account for the process of Syk de- 
phosphorylation in the model of Rcfs. ,33, 34_, or more 
generally, downstream influences on events triggered by 
scaffold-mediated nucleation. This feature of the ternary- 
complex model (e.g., omission of processes affecting Syk 
phosphorylation status, including dephosphorylation me- 
diated by phosphatases) explains the quantitative dis- 
crepancy between the two curves of Fig. [2JB). 
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B. Concentration-dependent scaffold functionality 

A bivalent scaffold, depending on its concentration, 
can recruit different pairs of ligands. In Fig.|3]A_, we show 
that two ligand pairs can be recruited to the same scaf- 
fold and to the same maximal extent at different scaf- 
fold concentrations. Co-localization of the first ligand 
pair is maximum at a scaffold concentration of 0.055 /iM, 
whereas co-localization of the second ligand pair is max- 
imum at 9.8 /iM. In general, such a scenario requires 
that the ligand pairs have different total concentrations 
and/or scaffold affinities (Fig. [3|3). The contour plot of 
Fig. [3J3 shows how C^'^'^, the concentration of ternary 
complex at Sq — S'p^*, depends on ligand concentration 
and affinity for the special case where Aq — Bq = Lq 
and Ka = Kj, = Kd. Thus, at different concentrations, 
a scaffold can produce quantitatively identical outputs, 
as measured by co-localization of its binding partners, 
but qualitatively distinct outputs, because the binding 
partners differ. Any scaffold is potentially multifunc- 
tional, as the concentration of a scaffold can be adjusted 
through regulation of gene expression. In addition, the 
effective concentration of a scaffold can be adjusted with- 
out changes in gene expression through post-translational 
modifications, e.g., phosphorylation-dependent binding 
sites on a scaffold can be turned on and off by kinases and 
phosphatases, effectively changing the amount of scaffold 
available for interaction with ligands. 



FIG. 3. A scaffold is potentially multifunctional. (A) Ac- 
cording to the ternary-complex model, a scaffold can recruit 
different ligand pairs at different scaffold concentrations. Pa- 
rameters used in calculations are as follows: Ao = Bo = 0.05 
/iM and iCa = -K't = 5 uM (solid line); Ao = Bo = 1.0 /iM 
and Ka = Kb — 8.8 /xM (broken line). Competition of ligand 
pairs was not considered explicitly in our calculations because 
competition is insignificant for this particular example. (B) 
A contour plot showing logiQ(C"i,'"'/l /iM) as a function of 
ligand concentration Lq (/iM) and equilibrium dissociation 
constant Kd (mI^)- We assume that Ka = Kb = Kd and 
Ao — Bo — Lq. Along any contour lin e, C *™"" is constant 



while Sg^' varies in accordance with Eq. ( 23 1 



C. Dependence of complex formation on binding 
constants 

Figure |4] illustrates how the maximum concentration 
of ternary complex, C™j^^, and the scaffold concentration 
at which Cab = C™^^ ('S'o^*) depend on the binding con- 
stants that characterize ligand-scaffold interactions. In 
Fig. [i]^, we can see that S^^^ increases as the equilib- 
rium dissociation constant [Ka = Kj, = Kjj for the 
results shown in this panel) increases, or equivalently, 
as affinity decreases. In Fig. we can see how ^q''* 
changes as the equilibrium dissociation constant for scaf- 
fold interaction with ligand A (Ka) is varied while Kb is 
held constant at 5, 50, and 500 nM. The value of Sq^^ is 
least sensitive to changes in Ka and Kf, for small values of 
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FIG. 4. Effects of binding constants on optimal scaffold con- 
centration Sq'^^ and maximum concentration of ternary com- 
plex C^i^^. (A) The value of Cab is plotted as the function 
of for three different values of Kd- In these calculations, 
Ka and Ki, are each set equal to Kd and Aq = Bo = 0.2 
fj,M. (B) The value of ^g^* is shown as a function of Ka for 
three different values of Kt- Ligand concentrations are the 
same as for panel (A). (C) The value of C^^" is plotted as a 
function of the ratio of the equilibrium dissociation constants 
(Kb/Ka)- In this panel, Ka is held fixed at 5 nM and ligand 
concentrations are the same as for panel (A). 



Ka and Kb (i.e., for high hgand-scafFold affinities). Fig- 
ure |4p shows how C^'^^ depends on Ka and Ki,. In this 
example, C^^^ is affected by changes of Ka and K^, but 
large changes of Ka or Kh result in only modest changes 
of C^^'^. For example, a 100-fold increase in Ki, from 5 
nM {Kb/Ka = 1) to 0.5 (Kb/Ka = 100) results in 
a only an approximate four-fold decrease of C^^'*'' from 
0.15 ^lM to 0.04 ^M. As indicated by Eq. the sen- 
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FIG. 5. Dependence of scaffold-nucleated complex formation 
on ligand concentrations according to the ternary-complex 
model. (A) The amount of ternary complex Cab is plotted as 
a function of total scaffold concentration So for three different 
values of Lq. In these calculations, Aq and Bo are each set 
equal to Lq and Ka = Kb — 5 nM. (B) Cab vs. So for three 
different total concentrations of ligand B (Bo) with Aq held 
fixed at 0.5 //M. The equilibrium dissociation constants are 
the same as for panel (A). 



sitivities of C^'^'^ and Sq^^ to an increase in a binding 
constant is analogous to the sensitivities of these quanti- 
ties to an increase in the concentration of a competitive 
inhibitor or decoy receptor. 



D. Dependence of complex formation on ligand 
concentrations 



Figure [5] illustrates how the maximum concentration 
of ternary complex, C^j^^^, and the scaffold concentration 

at which Cab = C'iT'^ (^'q''*) depend on ligand concen- 
tration s. T he results shown in Fig. [5] complement those 

I I _ I I I I. -opt 



of Eq. ([23]) and Fig. [s] In Fig. [H]^, we can see that S", 



increases as both ligand concentrations (^0 — Bq — Lq) 
increase simultaneously. Increasing ligand concentrations 
also increases C™j^^. Thus, inhibition of ternary complex 
formation by excess scaffold can be reversed by increas- 
ing ligand concentrations. For example, 10 /iM scaffold 
inhibits formation of the ternary complex at Lq = 0.5 
/iM. In contrast, when Lq = 1.5 /xM, C™^^ is increased 
by nearly three fold and there is a significant amount 
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FIG. 6. The effect of cooperativity on scaffold-mediated nu- 
cleation of ternary complex: Cab (mM) is plotted as a func- 
tion of 5*0 (m^) for five different values of the cooperativity 
factor (j>. The solid line corresponds to independent bind- 
ing of ligands A and B to the scaffold = 1). The broken 
line and dotted line below the solid line correspond to neg- 
ativity cooperativity: (j) ~ 0-1 ^'iid (j) = 0.01, respectively. 
Similarly, the broken line and dotted line above the solid line 
correspond to positive cooperativity: (j) = 10 and (p = 100, 
respectively. Calculations are based on the following ligand 
concentrations and binding constants: Aq = Bq = Lq = 0.5 



fj,M and Ka = Kb — Kd = 5 nM. In accordance with Eq. ( 23 1 
= Lo + Kd= 0.505 ijM for aU curves. 



of ternary complex at a 10 /iM concentration of scaffold 
(Fig. [sj'V). In Fig. [5|3, we show how Cab depends on the 
total scaffold concentration Sq for three different values 
of Bq (0.5, 1 and 5 /iM). For all three curves, the value of 
is the same (0.5 /xM). As can be seen, when the two 
ligand concentrations are different, the value of Cab can 
be nearly maximal over a broad range of scaffold con- 
centrations. For example, at Bq ~ 5.0 /iM and Ao ~ 0.5 
/iM, Cab is nearly maximal for scaffold concentrations be- 
tween 1 and 5 /iM. Over this range, ligand B is in excess 
and ligand A is limiting (i.e, nearly all A is bound to 
scaffold and the amount of ternary complex is approxi- 
mately equal to the concentration of ligand A). At higher 
scaffold concentrations, free B is consumed and Cab be- 
comes less than maximal because of the excess scaffold 
effect [TSKin]. 



FIG. 7. The effects of scaffold valence n on nucleation of a 
complex containing n ligands (Li,...,L„). The concentra- 
tion of complex containing all n ligands, Cn, is plotted as a 
function of the total scaffold concentration, 5*0, for three dif- 
ferent valences: n = 2, 3 and 4. Calculations are based on 
Eqs. Q and ([5| and the following parameter values: Li^ = 1 
/iM and Ki = 0.05 /iM for all i. 



F. Higher-order complexes 

So far, we have considered ternary complexes nucle- 
ated by a bivalent scaffold. Here, we consider scaffold- 
nucleated complexes containing two or more ligands per 
scaffold. In other words, we consider the reaction scheme 
of Eq. [l] which is illustrated in Fig. [T]\. In Fig. [Tj we 
show how C2, C3 and C4 depend on the total scaffold 
concentration Sq, where C„ is the concentration of scaf- 
fold bound to n ligands. As can be seen, for a larger value 
of n, the maximum value of C„ is reduced (in a controlled 
comparison). In addition, the scaffold concentration at 
which Cn is maximum is reduced. These results are con- 
sistent with the intuition that it is more difficult to form 
a scaffold-nucleated complex that contains more ligands. 
The results of Fig. [7] are essentially the same as those 
obtained by Heinrich et al. [TJ through the analysis of 
a more mechanistically detailed model, which provides 
further support for the relevance of the minimalist mod- 
els considered in this study. The asymmetry seen in the 
curves of Fig. [7| is explained by the asymptotic limits of 
|(/„|, which is given by Eq. ([t]), as 5*0 approaches and 



G. Effects of a scaffold on an enzymatic reaction 



E. Cooperative binding 

Cooperative binding of ligands A and i? to a scaffold 
is illustrated in Fig. ^ In accordance with Eq. ( 23 ) , S^^*" 
is independent of the value of (j). However, for ^ > 1, Cab 
takes on a nearly maximal value over a broader range of 
scaffold concentrations. For < 1, the maximal value of 
Cab is attenuated. 



In accordance with Eq. ( 30 ) , an enzyme-scaffold com- 



plex can be regarded as an effective enzyme, one that has 
a different maximal rate of reaction and one that achieves 
the half maximal rate at a different substrate concentra- 
tion, which can be lower or higher than the correspond- 
ing concentration in the absence of scaffold. Figure [8] 
complements these insights, showing how a scaffold that 
co-localizes a substrate and enzyme can affect the initial 
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FIG. 8. The effect of a scaffold on an enzymatic reaction 
under tlie assumption of rapid equilibrium. The initial nor- 
malized rate of reaction is plotted as a function of (A) total 
scaffold concentration So (/iM) and (B) total substrate con- 
centration Aq (/iM). In panel (A), = 1 A^M. The dash-dot 
curve corresponds to V/{kcatBo) (Eq. (29l), the dotted curve 



corresponds to V'/{hkcatBo) (Eq. (SOfyjand the solid curve 
corresponds to V/{kcatBo) + V /(hkcatBo). In panel (B), the 
dash-dot curve corresponds to the case of a classical Michaelis- 
Menten reaction (no scaffold). The other curves correspond 
to different scaffold concentrations: 0.5 piM (dotted line), 5 
jjM (broken li ne) a nd 5 (solid line). Calculat ion s are 
based on Eqs. |lO|)-( 15 I with (/> = I and with Eqs. (141 and 



( 15 I each modified to include a term for the concentration of 
{AB), which is taken to be AfBf /Km- The following param- 
eter values were used in all calculations; Ka = Kb = 0.005 
fj,M, Km = 0.6 fiM, and Bq = 0.2 fiM. 



rate of enzyme-catalyzed conversion of the substrate to 
product. As can be seen in Fig.[8j\, a scaffold can signif- 
icantly accelerate the rate of reaction relative to the sit- 
uation where the reaction proceeds in the absence of the 
scaffold. However, scaffold-mediated acceleration occurs 
only within a certain range of scaffold concentrations. 
A scaffold concentration outside this range, in the ex- 
cess scaffold regime, can have the effect of inhibiting the 
enzyme-catalyzed reaction. In Fig. |8j3, the rate of reac- 
tion is shown as a function of substrate concentration for 
different concentrations of scaffold and for the case of no 



scaffold. As can be seen, the dose-response curve can be 
shifted by the involvement of a scaffold, and in addition, 
the slope of the curve can be altered, i.e., the sensitivity 
of reaction rate to a change in substrate concentration 
can be altered. For the case of Fig. |8j3, the effect of the 
scaffold is to increase the sensitivity of the initial reaction 
rate to a change in substrate concentration. 



IV. DISCUSSION 

Elucidating the design principles of cellular regulatory 
systems is recognized as an important goal of systems 
biology [55H57j . which is facilitated by the tools of syn- 
thetic biology [38l [39] . Understanding the design princi- 
ples of scaffold-mediated nucleation of a multicomponent 
complex, a motif in cell signaling |14| . is important for 
a number of reasons, including understanding how scaf- 
fold function is conserved across species, interpreting the 
functional consequences of protein copy number varia- 
tion in diseases, and engineering scaffolds to have desired 
properties [40j. 

Here, we have used minimalist models for which ana- 
lytical results can be obtained to reveal principles that 
govern the recruitment of ligands to a scaffold. We 
have shown that the simplest model considered here, 
the ternary-complex model, captures essential features 
predicted by more sophisticated and mechanistically de- 
tailed models, which can only be analyzed via simulation 
(Fig. [2|. For the the ternary-complex model, we derived 
a number of analytical expressions, which can serve as 
basic design equations, in that these equations can po- 
tentially be used to guide the manipulation of intrinsic 
scaffold properties or scaffold milieu to achieve desired ef- 
fects on the system-level properties of scaffold-dependent 
processes. Perhaps the most significant analytical result 
is Eq. ( 23 ) . Special cases of this equation have been re- 



ported in studies of multivalent ligand-receptor binding. 
For example. Mack et al. [53] obtained a similar equa- 
tion in a study of bivalent ligand-receptor binding (cf. 
Eq. (16) in Ref. 29^ and Eq. (23) here in this report). In 
studies specifically focused on scaffolds, earlier work has 
largely relied on simulation. The parameter-dependent 
nature of simulation results can lead to incomplete un- 
derstanding. For example, on the basis of simulation 
results, Levchenko et al. [ini HI] reported that scaffold 
function is insensitive to binding constants, which is def- 
initely the case in certain parameter regimes. However, 
as can be seen from Eq. (23) and as illustrated in Figs.js] 
and [4| binding constants can have a significant impact 
on scaffold function. 

For the ternary-complex model, in our presentation of 
Eqs. (10)-(23), we have shown that the maximum con- 
centration of the ternary complex, C™j^^, and the scaf- 
fold concentration that yields this maximum concentra- 
tion, 5g^*, depend on ligand-scaffold binding constants 



and ligand concentrations (Eq. (23)). In addtion, C™j^^^ 



but not S'g'^', depends on the cooperativity of ligand A 
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and ligand B binding to the scaffold. Levchenko et al. 
[I6] reported that cooperativity shifts the optimal scaf- 
fold concentration (see Fig. 4B in Ref. [H]). This finding 
is inconsistent with the equilibrium binding properties 
of a scaffold (Eq. [23| . The shift of dose-response curves 
seen in the study of Levchenko et al. [16] is perhaps 
attributable to kinetic binding properties of a scaffold 
in the specific context of a MAP kinase cascade. The 
above mentioned analytical results are complemented by 
the simulation results of Figs. [3}|6] In our presentation 
of Eqs. we have shown that scaffold-mediated 

nucleation depends asymmetrically on total scaffold con- 
centration, which is illustrated in Fig. [7] This finding is 
consistent with the simulation results of Heinrich et al. 
[T7j (see Fig. 8 in Ref. [TTj). We also considered compet- 
itive ligands and decoy receptors, the presence of which 
can be mimicked by an increase in the equilibrium bind- 



ing constants for ligand-scaffold interactions (Eq. (27l). 
Finally, we found that a scaffold can reprogram the cat- 
alytic efhciency of an enzyme (Eq. ( 30 ) and Fig. [s]) . 

The results of Fig. [Sj'V suggest that scaffolds could have 
multiple context-dependent functions. A scaffold, de- 
pending on its concentration, can nucleate distinct sets of 
ligands. Thus, in a given cellular context, only a subset 
of a scaffold's known binding partners may be physio- 
logically relevant. In addition, in cases where the (ef- 
fective) scaffold concentration depends on the strength 
of a signal, a scaffold may enable a cell to respond dif- 
ferentially to signal strength. An example where effec- 
tive scaffold concentration depends on signal strength 
via post-translational modifications and other fast events 
(and not gene expression) is provided by the model of 
Goldstein et al. [53] and Faeder et al. [M] for signaling 
by the high-affinity receptor for IgE antibody (FceRI) . In 
this model, phosphorylated receptors in ligand-induced 
receptor dimers serve to co-localize two copies of the pro- 
tein tyrosine kinase Syk, which is necessary for Syk au- 
tophosphorylation and downstream Syk-dependent sig- 
naling events. Furthermore, the number of receptor 
dimers depends on the concentration of (multivalent) 
ligand that induces receptor aggregation. Thus, given 
the results of Fig. [3]^, receptor dimers in this system, 
which have a scaffold-like function, could nucleate differ- 
ent signaling complexes at different ligand doses [TJ |32] , 
with these complexes potentially eliciting different cellu- 
lar responses. We are not aware of a signaling system 
where scaffold concentration controls the formation of 
complexes that lead to distinct cellular responses, but 
given the number of binding partners of a typical signal- 
ing protein, we expect that scaffold multifunctionality, 
such as that illustrated in Fig. [3]A_, could be found. 

Results such as those shown in Fig. |4|C can be used 
to engineer a scaffold with non-native ligand-scaffold in- 
teractions to obtain a desired scaffold-dependent cellular 
response to a signal, assuming ternary complex formation 
correlates with the cellular response of interest. Park et 
al. [5] engineered the yeast scaffold protein Ste5 to inter- 
act with one of its binding partners, Ste7 or Stell, via 



a non-native protein-protein interaction. In these exper- 
iments, the Ste7 or Stell binding site on Ste5 was dis- 
rupted by a mutation of Ste5 (V763A/S861P or I504T) 
and Ste5 was fused to the PDZ domain of syntrophin. In 
addition, Ste7 or Stell was fused to the PDZ domain of 
neuronal nitric oxide synthase (nNOS). The nNOS and 
syntrophin PDZ domains heterodimerize with a Kd of 
0.6 [33]. In contrast, the native interaction between 
Ste5 and Ste7 or Stell is characterized by a Kd of about 
0.1 fiM [44l|45]. Using measured concentrations of Ste5 
(35 nM), Stell (39 nM) and Ste7 (68 nM) reported by 
Maeder et al. [33]) we can calculate the effect of increas- 
ing Ku from 0.1 to 0.6 ^M. We find that this change 
is expected to cause a 15-fold decrease in the amount of 
the ternary complex composed of Ste5, Ste7 and Stell. 
This prediction is consistent with the decrease in phos- 
phorylation of Fus3 observed by Park et al. 5 for the 
engineered Ste5-Ste7 or Ste5-Stell interaction. 

We have focused on equilibrium properties of scaf- 
folds that nucleate ternary and higher-order protein com- 
plexes. We also considered the effect of a scaffold on 
the rate of an enzyme-catalyzed reaction under a rapid 
equilibrium assumption. The results that we have found 
should be useful for understanding cellular responses 
that depend on the concentration and composition of 
a scaffold-nucleated complex, assuming that these re- 
sponses are insensitive to the kinetics of complex for- 
mation. As Locasale et al. [201 121] and Thalhauser 
and Komarova :22| have shown, the kinetics of scaffold- 
mediated complex formation can have physiological rele- 
vance. Thus, an important caveat of the results presented 
here is the assumption of equilibrium or pseudo equilib- 
rium. Another caveat is the assumption of monovalent 
ligands. If the ligands that interact with a multivalent 
scaffold are themselves multivalent, then complex dose- 
response behavior can arise, such that complex formation 
is maximal at multiple scaffold concentrations [T] [31] . 

An interesting future application of the results pre- 
sented here might be interpretation of comparative pro- 
teomics data. It has been observed that protein copy 
number tends to be conserved across species ^46ri48j . Be- 
cause the function of a scaffold is affected by the relative 
abundances of the scaffold and its ligands, conservation 
of relative abundances of the proteins in a "scaffold mo- 
tif" would suggest conservation of function of the motif, 
whereas the opposite would suggest the acquisition of a 
new function or compensatory modifications of scaffold 
properties and/or milieu, which could perhaps be more 
easily understood or predicted in light of this report. 
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